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Abstract 

In this paper, I present a technique to simplify the tensorial re- 
duction of one-loop integrals with arbitrary internal masses, but at 
least two massless external legs. By applying the method to rank I 
tensor integrals, one ends up with at most rank 1 tensor functions 
with the initial number of denominators, plus tensor integrals with 
less denominators and rank < L To illustrate the algorithm, I explic- 
itly compute diagrams contributing to processes of physical interest 
and show how the usual numerical instabilities due to the appearance 
of Gram determinants can be controlled. 
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1 Introduction 



Present and future experiments in high energy physics require an increasing 
understanding of multi-particle final states. Idealized calculations, in which 
heavy particles are treated as being on-shell, are often not accurate enough to 
match the experimental needs. In fact, the produced heavy particles always 
decay, giving a multi-body final state as an observable signal. 
At first, investigations of multi-particle processes can be performed at the tree 
level, but, in a second stage, also radiative corrections have to be taken into 
account. At LEP2, for instance, a complete electroweak one-loop calculation 
to off-shell four-fermion production beyond the leading log approximation is 
still missing, and this reflects on the accuracy in the W mass measurement. 
Besides, next-to-leading order QCD calculations of multi-jet final states are 
becoming unavoidable to control the background in the search of new sig- 
nals. Therefore, it is clear that multi-leg loop calculations are playing a role of 
increasing importance and that any progress achieved in this field is welcome. 

One of the major problems in loop calculations is the reduction (in n 
dimensions) of rank / tensor m-points integrals as follows 



qt^i . . . qi^i 



Di = {q + Si^i)^ - mj 

i 

■^i = H ko = 0. (1) 

j=0 

The standard reduction method ^ consists in decomposing the tensorial 
structure in terms of momenta kj and metric tensor. For example, a rank two 
5-point one- loop integral can be expressed in terms of 11 linearly independent 
scalar form factors 



= E2ik'^k^,+E22kl^k', + --- + E,ng'"'. (2) 

Analogously, the rank three 5-point function J^"^'' gives 24 terms. 

The freedom in choosing the basis for the tensorial decomposition can be 

used to define suitable combinations vj of kj such that Vi ■ kj = 6ij 0. In 
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any case, in order to express the coefficients of the decomposition in terms 
of scalar loop functions, one has to solve algebraical systems obtained by 
multiplying both sides of eq. by kj (or linear combinations of them) and 
the metric tensor. This causes the appearance of Gram determinants 

A = det{ki ■ kj) 

into the denominator, that can give rise to numerical instabilities. In fact, 
they vanish in coUinear regions of the phase space where the cross section 
is well defined, so that large numerical cancellations are expected in the 
numerator in order to keep everything finite. 

Different approaches also exist in which the form factors are identified with 
integrals in a different number of dimensions [§] or directly formulated in the 
Feynman parameter space In both cases. Gram determinants still appear 
at some stage in the reduction. 

Recently, a new method has been proposed by Campbell, Glover and Miller 
in which the coefficients of the reduction are built up by combining the 
scalar integrals in groupings that are well behaved in the limit A ^ 0. This 
somehow solves the problem of the Gram determinants, but to actually con- 
struct the well behaved groupings can be cumbersome in general - non QCD 
like - cases, especially when the rank of the appearing tensor integrals is 
large. Furthermore, the number of the basic functions increases, and addi- 
tional scalar integrals must be evaluated in n = 6 + e or higher dimensions. 

In this paper, I propose a method to simplify Feynman diagrams, in which 
no evaluation of new functions is required, apart from the usual scalar inte- 
grals (new functions may appear when n = 4 + e, but their computation is 
trivial). Although, in general, the appearance in the denominator of quanti- 
ties that can vanish in some corner of the phase space cannot be completely 
avoided, a better control of those singularities is possible with respect to the 
traditional methods of refs. |^ . The procedure is somehow complementary 
to the technique of Campbell, Glover and Miller, and can be used to simplify 
the problem before applying their method. 

A trivial example serves to illustrate the basic idea. Consider the tensorial 
reduction of the following quantity 



One can certainly use eq. (0), but the equation 

(^^ = = Di + ml 
immediately gives the desired answer in terms of just two scalar functions 

Therefore, by using in a diagram the algebra of the 7 matrices to reconstruct 
denominators, one gets a simpler tensorial structure. 

This technique of reconstructing denominators has been already used in the 
literature for the evaluation of specific integrals 0, 0]- In this paper, I show 
that a reduction based on this procedure can always be worked out for generic 
one-loop integrals with arbitrary internal masses, but at least two massless 
external legs. 

In the next section, I present the method, with the help two 4-dimensional 
examples. In section 3 and appendix A, I describe the extension to n = 4 + e 
dimensions and give explicit results for cases of physical interest. 



2 The Method 

After the decay of the intermediate heavy particles, a generic m-point loop 
diagram can be written as follows (see fig. 1-a) 

where pi...r are the r external momenta of the diagram, ki...m~i the momenta 
in the loop functions - as defined in eq. (P - and Tr*^*-* traces over 7 matri- 
ces, which contain ^'s. In most practical cases, the r external momenta are 
massless, in fact they can only be photons, gluons or light fermions. The 
appearance of traces in eq. (|^) is a general feature. In fact, if the external 
particles of the diagram are fermions, traces naturally appear, while the po- 
larization vectors of photons and gluons can always be written in terms of 
spinor strings using the following representation 0| 

e'^(X,k) = ^^^^ ^ b = arbitrary massless vector. (6) 
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In standard calculations, each term in eq. (^) is rewritten as follows 




and the tensors 



/ 



(Tq 



q^...qP ... 

Dl---Dm 



(8) 



decomposed using a basis of independent four-vectors, as described in the 
introduction. Instead, I propose to use 7-algebra and spinor manipulations 
to move all ^'s close to each other, and subsequently apply = q^ in order 
to rewrite each trace as 



where Tr^ contain at most one ^, Trf do not contain ^ and the coefficients 
Ak and Bi are functions of the scalar products between q and the external 
momenta of the diagram pi...r- Since a single ^ appears in the right side of 
6q- (H), just a rank 1 tensor decomposition is formally left. Anyway, there 
is no gain in applying the described procedure if the coefficients and Bi 
are generic. In fact, if scalar products of the kind {q ■ pj) appear and pj is 
not equal to one of the momenta kj in the loop functions, the same kind 
of tensorial decomposition necessary for computing eq. (|^) is hidden in eq. 
(P). On the other hand, if only powers of (g ■ kj) and appear, it is always 
possible to express at least one of them in terms of denominators. In this 
case, starting from m-point rank / tensor integrals, the algorithm gives at 
most rank 1 m-point functions, plus n-point rank p tensor integrals with 
n < m and p < I. 

Naturally the question arises whether the coefficients A/^ and Bi can be made 
to have such a nice feature. To answer that question, while keeping every- 
thing transparent, I will assume, for the rest of this section, four-dimensional 
space time, namely ultraviolet finite integrals, with four-dimensional regu- 
larization of infrared and collinear divergences. The extension to r?, = 4 + e 
dimensions will be considered in the next section. In the following, I intro- 
duce all notations and necessary formulae. 



Tr(^) [^...f..]=Y^A, Trf • •] + E Tr? [■ ■ ■] 



(9) 



k I 
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Basic objects are strings of 7 matrices between massless spinors. I will use 
the notation 

^+ (*) 1^3 ^ 1^5 ■ • ■ ^ (j) = {3 4 g 5 ■ ■ ■ 6}+T . (10) 
The Weyl spinors are defined as 

u±ij) = n± u{j) u±{j) = uii^ 

v±ij) = n^vij) v±ij)=vU±, (11) 

where j denotes the spinor momentum and Il± = |(1 ± 75) are chirality 
projectors obeying 

n+n_ = o (n±)2 = n±. (12) 

Note that, using u±{j) = %(i), one can always rewrite any spinorial string 
with V spinors on the left and u spinors on the right, so that the notation in 
eq. (|l^) is completely general. 

The completeness relations for massless spinors read 

u±{j)u±{j) = ^jU^ 

v±ij)v±ij) = ^Hi, (13) 
from which one can easily get similar relations for generic momenta i and j 



u^{i)v-{3) = -t^fVy"-' (14) 

where b is an arbitrary four- vector and 

<ij>* = v_{i)u-{j) = u+{i)v+{j) 

<ij> = -v+{i)u+{j) = -u^{i)v_{j) 



(15) 
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are scalar products in the spinor space p|, ITO 



The following identities hold for spinor strings |11] 



{r} 



i 3 



{r}. . , 



(16) 



r being the string obtained from F by reversing the order of the 7 matrices. 
The following Kahane-Chisholm identities [0 are also very useful 

TaFoT" = -2 f „ 7,re7" = Tr[T,] - 75 TrYi^T,] , (17) 

where Fq and Fe stand for strings with an odd and even number of 7 matrices, 



respectively. Finally, one can prove 



if 6^ = 



(18) 



Armed with the previous formulae, I will show that, in four dimensions, the 
coefficients Ak and Bi in eq. can be made to depend only on powers of 
and (q-kj) on the condition that at least two momenta in the set ki...m-i (say 
ki and k2) are massless, namely coincide with some of the external massless 
momenta of the diagram, as in fig. 1-b and 1-c. 
In fact, the following identity 






Pr Pb Pr P5 Pr PS 

a b c 

Figure 1: Generic one-loop multi-leg process (a). The blob stands for the sum of 
all possible one-loop diagrams. Two classes of diagrams are also shown, for which 
the reduction method is applicable (b and c). 
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allows to "extract" q from the traces 



Tr[T^iT2] = — —{2{q-k2)Tr[T4,T2]+2{q-kr)Tr[T42T2] 

2 [ki ■ k2) 

- Tr[TMh^2] - Tr[TiMh^2]} , (20) 

where Fi 2 represent generic strings of 7 matrices. In the first two terms, the 
q dependence is already extracted out. The other two traces can be broken 
with the help of eq. ([T3|) 

Tr[r,h02r2] = Tr[ri^i#2(n+ + n.)r2] 

= {r2ri}2+T + {q};1, {r2ri}-t • (21) 

By iteratively applying the above procedure together with the identities (de- 
rived again by means of eq. (^) and eq. (|l4|)) 

{qKUiK^i = ^{q-h){q-k2)-2q\h-k2) 

WrtWrt = i4Ty^^[2('?-^i)M2^T-2(9-6){i2g}+T 

+ 2(g-A;2){16g}+7-gMl62}+7] (22) 

it is easy to see that the resulting coefficients and Bi only depend on 
powers of g^, (g ■ ki), {q ■ k2) and {q ■ b). The vector b is arbitrary (and not 
necessarily massless), but different from ki and k2. By identifying it with a 
third momentum in the set ki...m-i Q, Ak and Bi have the suitable structure 
to get the desired simplifications in the tensorial decomposition. 

In the computation of physical processes, the one-loop diagrams with 
highest m always lie in the class of the corrections connecting two external 
legs (see again fig. 1-b and 1-c). That guarantees at least two massless 
momenta in the loop functions and, therefore, that the described reduction 
method is apphcable at least for the most complicated diagrams appearing 
in the calculation. 

Finally, one should notice that structures that can vanish appear in the 
denominator of eqs. (|20|) and (|2^). However, as we will see, often in practical 

^This identification is possible only for m-point loop functions with to > 3. However, 
in the next section, I will show that the method still works when m < 3. 
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calculations there is no need to extract the q dependence as described, since 
the reconstruction of denominators naturally takes place without introducing 
too many factors in the denominators. 

In the following two examples, I show the method at work in four dimensions. 



Example 1 

I shall compute, in the renormalizable gauge, the 6-point diagram of fig. 2, 

relevant for studying electroweak corrections at LEP2. 

The spinorial structure in the numerator of the diagram reads 

= -{liB(3}tAPMtA^^Cv}t2 , (23) 

from which the scattering amplitude M. can be obtained after division by 
the denominators and integration over d^q 



M = j (fq 



A 



Di D2 D-i D^D5 De 



pa = q pb = q + P4 + P5 
Pc = q-P2-p& 

D, = q' 

D2 = {q + p,f-Ml 

D-i = {q + P^ + P^y 
Di = {q + P5+Pi+ pzf - 

D5 = {q+P5+P4+P3+ Plf 
D6 = {q+P5+Pi+P3+Pl+ P2? - 

Figure 2: LEP2 6-point diagram. All momenta are incoming. V can he either a 
massive gauge boson or a photon. In the latter case, an appropriate 
four- dimensional infrared regularization is understood. 

Using eq. (0) one can rewrite 

A = -^{^B(34b5(^Au6blfiCu}^^2 
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-— {166z/A4655Cz/}+2 



4 

N 

N = v+{4:)pu_{5) v+{6)^u4l), (24) 



where I have apphed twice the Chisholm identity in eq. (0) to sum over P 
and /i. By introducing a new summed four-dimensional index a, the products 
u_{i)v^{j) can be reconstructed 

A = 2 v43h^u46) v+il)rYhu-i4^) v+i5)McluU-i'^) (25) 

and, after Fierz reordering the last two strings by means of the second of eqs. 
(O), one gets 



A = 4i;+(3)7«iz_(6) t;+(l)7"«_(2) ^+(5)^|^c|^^ii_(4). (26) 

Since 

i>Bi>ci>A = [D5 + 1^131^26]^ - Di 

= h+i>3 (27) 

denominators naturally appear in A^ so that the scattering amplitude is 
immediately expressed in terms of one scalar 5-point function, one rank one 
5-point function and one 6-point function or rank 1. 

Example 2 

A different contribution comes from the "crossed" 6-point diagram in fig. 3. 




Figure 3: "Crossed 



q Pb = q + P3+P6 

q-pi-pA 

q' 

{q + Ps + PeY 

{q + P3+P6+P5? - 

{q+P3+P6+P5+P2Y 
{q+P3+P6+P5+P2+ Pi? - 

EP2 6-point diagram. 
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The numerator of the amphtude reads 

A = {z/5/?}+- . 

Exactly as before, one easily performs the sums over ^u, p and (3 



(28) 



A 



-8 {B}t2 {C]t, = 8 {g}^6 {6}^2 {1} 



I 5 4 



{9}^6 {6}^2 {^}^; + {9}^6 {^}3 2 {1}^4 - {^}^6 {^Is 2 {?} 



5 4 



(29) 



The first term in eq. ( |29| ) only contains one g and cannot be further manip- 
ulated, while the second and the third term can be reduced with the help of 
eq. (0). For example 



4 5 



< 64 > 



{g64g}+ 



2 {q ■ P,) - 2 (? ■ Pa) mtt + <f mtX ■ (30) 



< 64 > 

Similarly, one gets for the last term 

{9)3^2 = W^6{9}2:W^4 

1 



{g62g35g}+ 



1 



< 62 >< 35 

< 62 X 35 >* ■ " P'^) {25g}^ 4 - 4 (g ■ Pe) (g ■ Ps) {23g}+- 
+ 2 (g ■ pe) g2 {235}+4 - 4 (g ■ p^) {q ■ Ps) {65g}+; 

+ 4 (g ■ P2) (g • Ps) {63g}n - 2 (g ■ pa) {635}+^ + q' {6235g}+4 ] . (31) 

Therefore, by reconstructing denominators from the dot products, one easily 
see that the most complicated objects appearing in the scattering amplitude 
are rank one 6-point functions and rank two 5-point functions, as expected. 



Some comments are in order. What I have done in the previous examples 
is equivalent to Fierz reordering the spinorial strings. This already gives the 
correct structure to apply the reduction algorithm, without using eq. (pJ|). 
The fact that some diagrams are simpler to compute than others, depends 
on what one gets after Fierz reordering. In fact, after elimination of all sum 
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indexes, one can end up with momenta containing q (that I denote with 
capital letters) either close to each other or "sandwiched" between spinors. 
Namely 



{ABC---y^ 



or 

{A}tj{B}t-, • 

In the first case, the final result is simpler because the equation = can 
be immediately applied, while, in the second case, one first needs to use eqs. 
(|14D and some 7- algebra. 

In example 1 there are no singular structures, while, in example 2, a single 
quantity, that can vanish in the collinear limit, appears in the denominator 
(see eqs. (|30| ) and (PI)). This roughly corresponds to a situation with a single 



Gram determinant in the denominator, that is known to be numerically quite 
stable p (actually, since < ij >< ij >*= 2{pi ■ pj) eqs. (|30D and (^) are 
even more stable). In addition, even in this second case, one can in principle 
get rid of all singularities by changing the reduction when the collinear limit 
is approached. For example, a reduction equivalent to eq. (PT|), but with a 
different singular structure, can be obtained by replacing 

1^3 2^4 3^5 

4^6 5^1 6^2. (32) 

In conclusion, one has a better control on the singularities, with respect to 
traditional reduction methods, in which one would have immediately to start 
with rank three 6-point functions and a higher number of Gram determinants. 



3 The n-dimensional case 

In this section, I show how the method can be applied in n dimensions. 
The basic idea is simple. Starting from the traces in eq. (|^) one breaks each 



n-dimensional vector into four- dimensional and e-dimensional part 0, |T4 
(e = n — 4). This allows to factorize any e dependence out of four-dimensional 
objects, in which all properties of the spinor calculus can still be used in order 
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to work out the reduction, exactly as shown in the previous section. This 
procedure is equivalent to assume a four-dimensional-helicity regularization 
scheme, in which only the loop momentum is continued in 4 + e dimensions, 
while helicity states are kept four-dimensional ||15|| . 

Before illustrating the procedure with explicit examples, I introduce nota- 
tions and relevant formulae. 

I underline quantities that live in the n-dimensional space, while a tilde 
denotes e-dimensional objects. A generic n-dimensional vector v can then be 
written as 

v = V + V . (33) 

It is always possible to chose orthogonal spaces for 4 and e-vectors, such 
that Vi ■ Vj = 0. Besides, e-dimensional 7 matrices freely anticommute with 
four-dimensional ones 

{7m, 7.} =0. (34) 

As for 75, a mathematical correct prescription would be keeping its definition 
also in n — 4 dimensions, that implies 

[75,%] =0. (35) 

While this causes no trouble for diagrams with only infrared or collinear 
divergences, extra finite terms, that violate the Ward identities, appear when 
ultraviolet divergences are present. At the one-loop level, one may either 
adopt eq. ( p5D and restore the Ward identities, by introducing appropriate 
finite counter terms in the Lagrangian, or impose, also for e-dimensional 
matrices 

{75,%} = 0, (36) 
plus extra care in manipulating traces [|TB1. In the following, I always assume 



eq. (pq). To preserve consistency, it is enough anticommuting 75 to the end 
of a trace before performing any 7-algebra [|17|. 

As a first example, I compute the simple vertex diagram in fig. 4. 
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w 




pa = q + pi 
pb_ = q-p2 

D2 = {q + Pif 
D, = {q-p,f 



Figure 4: Vertex diagram in n dimensions. 



Splitting up the n-dimensional traces as in ref. |T3], the numerator of the 
amphtude can be written as follows 



1 2 



e {A^iBY,- - + e~q' {^j^" , (37) 



where all e-dimensional 7-algebra has been explicitly worked out in order to 
get the desired factorization of the e-dimensional objects. Summing over /? 
gives 



A^ 



-2 2 - e {qm}r2 + r (2 + e) {/x}; 



2 ' 



(38) 



where the Dirac equation has also been used. 

In the first term of eq. (pSj) , A and B can be extracted with the help of eq. 

m 



1 



4 {pi ■ 



^{PA-P2) (pb-Pi) 



-2 {pB ■ Pi) {2^2gl}+2 - 2 (pa ■ P2) {2qlfil}t2 + {2qlf^2ql}l 



I +- 



(39) 



One easily convinces oneself that the last three terms of eq. ( p9D do not 
contribute to the amplitude. In fact, performing a four-dimensional tensorial 
decomposition, each coefficient gives zero. Therefore, one gets the following 
result in terms of denominators 



D1-D2 

Pl-P2 



(40) 
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Analogously 

and the final result reads 

= {/i}+- 



{2 + e) (f + ^^^{q-p.) 



(41) 



Pi ■ P2 

4 (pi ■p2) + 2 (D3 + D2-2 Di)] . 



(42) 



In the previous equation, the desired reconstruction of denominators has 
been achieved, but a term proportional to is also left, that requires the 
evaluation of the extra scalar integral 



(Tq 



By decomposing the integration as follows [|18 
J (Tq = J d^q d^jj with 

one easily gets 



TT 



(43) 



(44) 



(45) 



that is all one needs. 

The result in eq. (^2]) can be checked by using standard reduction techniques. 



As a second example, I compute the box diagram in fig. 5, contributing to 



the one- loop QCD corrections to 4 fermion production at LEP2 |]T9|. 




pa = q + pi PB_ = q-P2 
Pc = q + Pi + P3 

D2 = {q + Pi)' 

D3 = (q + Pi+Ps)^ 

D, = {q-P2f 



Figure 5: QCD one-loop box diagram. 
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The spinorial structure reads 

A""^ = - A a C n B 

= 2 {BuCaA}^^ + e {qaCnq}^^ 

- 2f\{fiaA}^^ + {Bfia}^:2+{iiCa}^ 



Decomposing the integration as in eq. (@), one gets 



rq 



(f q^ 



+ O{e) 



In addition 



/ 



d^'q 



{qaCfiq}^^ 



(46) 



(47) 



(48) 



does not contain any pole. In fact, infrared divergences are absent by power 
counting and, in the kinematical regions where coUinear divergences take 
place - namely when q (x pi or q oc p2 - the numerator of eq. (|48|) vanishes 
thanks to the Dirac equation. Therefore, in the limit e ^ 0, only the first 
term in eq. p6|) survives 



A 



afj, 



2 {BfiCaA}^ 



(49) 



Notice that, since the diagram is free from ultraviolet divergences, this result 
coincides with what one would get by assuming eq. (|35|), as expected. 
Terms containing up to three powers of q appear in eq. (|49|) . One could 
in principle extract the q dependence by means of eq. (^), but this would 
cause the appearance of (pi ■ P2) into the denominator. This can be avoided 
as follows. Starting, for example, from 



2 {g^gal}+2 
one gets, after two anticommutations 



(50) 



2 {qfiqal}^ 



4{a}+2 qYi-2{q-p,)q^ + Ap'^ {qfiq}^' (51) 
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and, rewriting the last term trace 

- 1 If 

{<lM}t'2 = ^Tr[g231g//n+] = — [2 (g ■ pa) {g/^3}+7 

- 2 (q- ps) {qfi2}^- + 2 (g • pi) {fiq3}+', - q^ {32fi}+- 
N = {3}+7, (52) 

from which the reconstruction of denominators can be performed as usual. 
In eq. (|52D, appears into the denominator, instead of (pi ■p2)- N does not 
vanish in the collinear limit pi (x p2, but gives zero when p^ oc pi or p^ oc p2- 
However, since ps refers to a massive particle, those configurations practically 
never occur in numerical programs, so that the described reduction gives 
numerically stable expressions. 

All terms in eq. (^Up can be reduced as shown. The complete result is pre- 
sented in appendix A. Notice that the reduction can be worked out exactly 
in the same way even when one of the internal fermion lines of the diagram 
is massive. Only denominators change. For example, with the replacement 
(pi + Ps)^ — >■ (pi + Ps)^ — mf, the result in appendix A can also be used to 
compute the box diagram that one encounters in evaluating QCD corrections 
to off-shell single top production at hadron colliders 



As a last example, in order to introduce the algorithm in the presence of 
external massless vector particles, I show the reduction of the QCD four- 
point function in fig. 6, appearing in the calculation of the W decay into 3 
jets. 

D, = q' 

^ /\^^ D2 = {q + Pi? 

^ ^ W^3, a ^3 = (^_p^_p3)2 



Figure 6: Box diagram inW^S jets. 
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As before, all terms in e dimensions do not contribute, when the limit e — 
is taken, so that only four-dimensional objects appear in the numerator of 
the diagram 



yfiaS 



/3 



(53) 



The polarization vector of the gluon can be expressed in terms of spinors by 
means of eq. (Bf). For example, by choosing h = pi when A = — 1, one gets 



I J3 1 



e„(-l,p3) = 

= < 13 >* . 
Therefore, eq. (|53D can be rewritten in terms of single traces 



(54) 



A^ 



{pi q fi {q - P2 - P3) iq + P3)} 



3 2 



+ {(?-P2)^(? + Pi) (g-2p3)pi} 



3 2 



2 {qpi {q-p2- Ps) t^iq + Pi)} 



3 2 



(55) 



and the reconstruction of denominators performed in the usual way. For 
completeness, I list the result in appendix A. 



4 Conclusions 

I have presented a method to simplify the tensorial decomposition of loop 
diagrams in n dimensions. Rank / tensor integrals, with at least two massless 
momenta, can be reduced to at most rank one tensor functions with the initial 
number of denominators, plus integrals with less denominators and rank < /. 
A better control of the numerical instabilities due to the appearance of Gram 
determinants is possible, with respect to traditional reduction techniques. 
For processes where all heavy particles decay, the class of the most compli- 
cated loop diagrams lies in the domain of validity of the proposed method, 
allowing strong simplifications. 

To illustrate the approach, I computed diagrams contributing to several pro- 
cesses of physical interest. 
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Appendix A 



In this appendix I present the complete reduction of the diagrams in fig. 5 
and 6. 



The amplitude of the diagram in fig. 5 is 



M 



(Tq 



Di D2 Ds 



{A.l) 



with A°''^ given in eq. (^91) . 

One gets, in terms of denominators and = {3}^^, the following expression 



A' 



{qfi{l + 3)al}+2 - {2fi{q + 1 + 3)al}+2 - {2/i(l + 3)ag}+ 



2 

+^ {(^2 - ^1) [{/ig3al}+2 - {9/^31a}+2 - {aq32fi} 



+- 
2 2 



+ {aiq + 1 + 3)/ig3}+2 + -^^^^^ {{a23qfi}+- + {a23/ig}+7 

- {/ig32a}+;)] + (Di - D^) [{g/i31«}+7 + {ag/i23}+; - {ga2/i3}+; 

+ {qi,{q + 1 + 3)a3}+7 - i^l±^ {alg/.3}+7 

+Di - 2)a23}+7 - {32/i(l - 2)a}+7 - {32a(g + 1 + 3)/i}+7 

(Pi+Pa)^ 



N 



({a231/i}+i-{al32^}+7 



+ (A3 - Do 



{ga2^2}|; 



- {g//21a}+; - {qii{q + 1 + 3)a2}+7 + ^Il±pl. {alg^2}2+; 



+ {Pi + Vzf {g/"21a}+7 - {(za2/i2}+; + {q^i{l + 3)a2}+; 
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N 



{alg//2}+^ 



(A2) 



The above equation has been numerically checked against the result in ref. 



1^, that was obtained using the standard reduction techniques of ref. 
As for the diagram in fig. 6 (with gluon polarization A = —1), one has 



AH-V, 



(A3) 



with 



A^(-l) = ^{{2D,-D,){lq^i}^- + 2D,{l^,q}^,- 

+ (D2~A{p,-ps)) {(g-2)/il}+2 
+2 (D, - D2) [{qfi{q + 1)}+- + {2qf,}^~, 

+2Di [{l/il}3+2 - {12/i}+2l +2 {2/ig31}+- 



+2 {gl(2 + 3)/xl}+2} 



= < 13 >* 



(AA) 
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